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PUBLICATION  REVIEW 


ABSTRACT 


This  report  summarizes  the  results  of  an  investigation 
which  was  carried  out  during  the  period  January  1,  1972  to 
December  31,  1972.  The  purpose  of  this  study  was  to  investi¬ 
gate  the  potential  advantage  which  may  accrue  to  HF  array 
systems  by  employing  adaptive  beamforming  techniques.  Data 
for  the  investigation  were  obtained  using  a  2600  km.  path 
between  a  CW  transmitter  located  at  Bearden,  Arkansas  and  a 
receiving  array  at  Los  Banos,  California.  These  data  were 

Originally  recorded  for  purposes  of  measuring  the  azimuthal 
resolution  of  the  ionosphere. 

The  data  consisted  of  digitized  recordings  taken  sequen¬ 
tially  from  eight  individual  whip  antennas  in  a  linear  array. 
Elemental  spacing  was  360  m.  providing  a  total  aperature  of 

about  2.5  km.  a  receiver  bandwidth  of  31.25  Hz  was  used  for 
each  of  the  eight  elements. 

For  the  purposes  of  this  research,  the  Bearden  signal  was 
presumed  to  be  an  undesired  interference.  The  array  element 
outputs  were  combined  through  tapped-delay-line  filters  for 
purposes  of  rejecting  this  interference  while  simultaneously 
forming  a  beam  in  the  direction  of  a  weak,  desired  signal.  It 
was  shown  that  by  adaptively  adjusting  the  delay-line  weights, 
the  Bearden  signal  could  be  reduced  by  more  than  25  db  over 
conventional  time  delay  and  sum  beamforming.  The  3  db  adaptive 
time  constant  was  measured  to  be  about  27,000  adaptations, 
corresponding  to  a  real  time  rate  of  less  than  100  msec  for 
high-speed,  special-purpose  digital  processors. 
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I.  SUMMARY 

This  report  contains  the  technical  results  obtained  under 
RADC  Contract  No.  F  30602-72-C-0 386  during  the  period  January  1, 
1972  to  December  31,  1972.  The  primary  results  are  summarized 
briefly  in  this  section.  Detailed  descriptions  are  presented 
in  the  sections  following.  The  objective  of  this  research  pro¬ 
gram  is  to  study  and  investigate  the  potential  advantages  which 
may  accrue  to  HF  antenna  array  systems  by  applying  adaptive 
signal  processing  techniques  similar  to  those  currently  used  in 
sonar  and  seismic  systems.  Data  for  this  study  are  obtained 
from  phase-coherent,  digitized  recordings  taken  using  high- 
frequency  radio  research  facility  at  Los  Banos,  California. 

The  recordings  consisted  of  data  samples  from  eight  individual 
whip  antennas  in  the  array  uniformly  spaced  at  360  meters. 

These  data  were  supplied  by  the  Ionospheric  Dynamics  Group  at 
Stanford  Research  Institute,  Menlo  Park,  California. 

On  March  15,  1972  two  digital  tape  recordings  were  re¬ 
ceived  from  SRI.  All  data  were  taken  over  the  Bearden, 

Arkansas,  to  Los  Banos  2600  km  path  under  conditions  of  high 
signal-to-noise  ratio  (25  db)  at  the  individual  whips.  Although 
these  data  were  originally  taken  to  investigate  the  resolution 
of  a  large-aperture  HF  array,  they  were  also  amenable  to  adapt¬ 
ive  processing  studies  in  that  the  signals  were  recorded  using 
eight  phase-coherent  receivers. 

Tt  was  decided  that  in  the  scope  of  the  present  research 
program,  the  Bearden  transmitter  would  be  treated  as  an  unwanted 
jamming  signal.  Adaptive  processing  was  then  used  to  extract 
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a  weak  signal  which  was  incident  on  the  array  from  a  slightly 

different  angle  and  at  a  different  frequency  than  was  the 

Bearden  signal.  Results  of  this  processing  were  then  compared 

with  those  of  a  conventional  array  processor  (time  delay  steer 

mg  with  amplitude  shading)  designed  to  receive  the  same  weak 
signal . 


The  reason  for  characterizing  the  Bearden  signal  as  un 


desired  interference  was  that  in  high  signal-to-noise  ratio 
environments,  optimal  array  processing  provides  approximately 
the  same  output  SNR  as  does  conventional  processing.  Under 
conditions  of  low  input  SNR,  however,  the  improvement  can  be  as 
high  as  40  db.  since  the  purpose  of  this  research  program  was 
in  estigate  the  potential  advantages  of  adaptive  processing, 
a  low  SNR  environment  was  selected. 

It  was  found  that  when  a  weak  signal  d1  displaced  in 
frequency  by  25  Hz.  and  in  angle  by  0.51  degrees  from  the 
Bearden  trinsmitter  was  selected  as  the  desired  signal,  adaptive 
processing  provided  an  improvement  of  29  db  over  conventional 
beamforming.  When  a  signal  d2  at  the  same  frequency,  but  dis¬ 
placed  in  angle  by  1.02  degrees,  was  used  as  the  desired,  the 
improvement  was  about  6  db.  These  results  are  summarized  below: 

COMPARISON  OF  PROCESSED  OUTPUT  SIGNAL-TO-NOISE  RATIOS 
PROCESSING  METHOD  OUTPUT 


SNR 

d_ 


individual  whip 
conventional  processor 
16  coefficient  adaptive 
40  coefficient  adaptive 
80  coefficient  adaptive 


-51.6  -51.6 
-43.3  -30.9 
-37.5 

-19.4  -25.0 
-14.4 


? 


The  adaptation  time  constant  for  these  studies  was  measured  as 
680  adaptations  for  the  40  weight  processor.  Existing  special 

purpose  digital  hardware  could  achievi  this  time  constant  in 
less  than  100  milliseconds. 

A  second  result  which  has  been  obtained  from  these  data  is 
the  design  of  a  simple  filtering  procedure  to  remove  equipment- 
generated  spurious  signals  from  the  raw  data.  This  procedure 
uses  the  calibration  signal  as  a  reference  and  has  been  shown 
to  reduce  spurious  levels  by  15  to  25  db. 

A  third  area  was  also  investigated.  The  adaptive  tapped- 
delay  lines  used  to  form  a  beam  output  have  previously  been  im¬ 
plemented  as  non-recursive  filters  -  that  is,  no  feedback  con¬ 
nections  are  employed.  It  is  well  known  that  feedback  delay 
lines  can  achieve  a  wider  range  of  transfer  functions,  for  a 
given  number  of  taps,  than  can  non-feedback  lines.  Feedback  has 
not  been  previously  implemented  into  adaptive  systems,  however, 
because  of  the  problem  of  maintaining  filter  stability.  During 
the  course  of  the  present  research,  a  procedure  was  devised  for 
adaptive  feedback  systems  which  is  guaranteed  to  produce  a 
Stable  filter.  This  result  is  new  and  represents  a  significant 
step  in  adaptive  technology.  Results  are  presented  illustrating 

this  procedure  for  the  case  of  a  single  input,  single  output 
filter . 
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II. 


INTRODUCTION 


A.  PURPOSE 

The  purpose  of  the  research  described  in  this  report  was  to 
investigate  the  potential  advantages  which  may  accrue  to  an  HF 
antenna  array  system  which  employs  adaptive  beamforming  techniques. 
The  direction  taken  in  this  study  was  to  compare  adaptive  array 
signal  processing  methods  similar  to  those  presently  used  in 
several  sonar,  seismic,  and  microwave  systems  with  conventional 
array  processing  methods  currently  employed  in  HF  systems.  Con¬ 
ventional  methods  generally  involve  the  use  of  time-delay  steer¬ 
ing  amplitude  shading,  such  as  Tchebeychev,  at  each  array  element. 
While  the  present  study  was  limited  to  a  comparison  of  single 
beam  systems,  the  results  obtained  are  also  applicable  to  mul¬ 
tiple  beam  comparisons. 

B.  MOTIVATION 

Recent  years  have  demonstrated  an  ever-increasing  number  of 
users  of  the  HF  spectrum,  and  a  resulting  increase  in  the  inter¬ 
ference  levels  experience  by  any  one  user.  Such  interference 
appears  as  unwanted  noise  in  the  users'  receiver  and  therefore 
degrades  his  siqnal-to-noise  ratio.  From  the  viewpoint  of  a 
communications  channel  this  degradation  is  reflected  directly  in 
terms  of  reduced  reliability  of  data  transmission  which,  in  turn, 
necessitates  a  reduced  information  rate  to  recover  reliability. 

In  order  to  improve  the  performance  of  HF  systems  operating 
in  a  greatly  increased  interference  environment,  users  have 
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resorted  to  three  basic  techniques: 

1.  increased  transmitter  power. 

2.  complex  modulation  and  siqnal  design. 

3.  increased  antenna  directivity 

The  first  of  these  is  clearly  the  least  desirable  since  an 
attempt  to  increase  s.ignal-to-noise  ratio  in  one  system  raises 
the  interference  level  in  all  other  systems.  The  second  tech¬ 
nique  has  become  increasingly  important  in  recent  years.  It 
involves  the  use  of  time  and  frequency  dispersed  signals  which 
employ  pulse  compression  techniques  at  the  receiver.  An 

example  of  this  approach  s  the  use  of  sweep-f requencv  signals 
in  ionospheric  sounders . 

The  third  technique  listed  above  is  the  one  most  closely 
related  to  the  present  study.  Previous  investigators  have 
reported  on  the  use  of  large  aperture  transmit  and  receive 
antennas  in  HF  systems.  The  work  described  in  Sweeney  for 

example,  showed  that  a  2.5  Km  receiving  aperture  antenna  can 
be  effectively  employed  in  the  HF  environment.  The  use  of 
narrow  beams  is  particularly  attractive  because  in  contrast  with 
increased  transmitter  power,  increased  antenna  directivity  in 
one  HF  system  will  have  the  effect  of  reducing  the  interference 
level  in  many  other  systems. 

The  research  program  described  in  this  report  was  undertaken 
to  investigate  a  fourth  and  equally  desirable  technique  for  inter¬ 
ference  elimination  —  namely,  the  use  of  automated  beamforming 
to  produce  directivity  nulls  in  the  direction  of  interfering 
sources  at  the  receive  antenna. 
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C.  BACKGROUND 

^he  research  effort  described  in  this  report  has  relied 
upon  the-  extensive  results  obtained  in  the  area  of  adaptive 
beaniforminq  in  the  fields  of  sonar,  seismoloqv  and  microwave 
systems.  A  complete  literature  survey  of  these  results  has 
been  completed  and  the  reader  is  referred  to  Reference  2  for 
detailed  information.  The  followino  paragraphs  present  a  brief 
summary  of  the  major  results  achieved  to  date. 

Since  the  work  of  Mermoz[3]  in  the  late  1950's,  a  signi¬ 
ficant  number  of  results  have  been  obtained  relating  to  the 
problem  of  improving  the  quality  of  signals  received  by  a  spa¬ 
tial  array  of  sensing  elements.  The  early  efforts  in  this 
development,  including  those  of  Mermoz  were  directed  toward  the 
processing  of  sonar  signals.  Shor[4},  Stocklin[5],  Vanderbulk [ 61  , 
Anderson^  ^ ,  and  Brynn ^  all  suggested  procedures  for  contain¬ 
ing  the  outputs  of  an  array  of  sensor  through  a  system  of  tapped- 
delay-line  filters  similar  to  that  shown  in  Figure  1.  This 
system  represented  a  marked  departure  from  conventional  array 
processing  which  consists  of  delayinq  the  sensor  outputs  in 
order  to  steer  the  array  in  a  prespecified  direction  and  then 
summing  the  delayed  outputs. 

Tapped-delay-line  processors  are  easily  constructed  and 
provide  a  wide  variety  of  spatial  and  temporal  processing 
flexibility.  For  example,  by  correct  choice  of  the  filter 
coefficients,  the  processor  can  null  out  several  directional 
interferences  while  simultaneously  presenting  a  preformed 
beam  in  the  direction  of  the  desired  signal.  Considerable 
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effort  has  been  devoted  to  the  development  of  optimality  cri¬ 
teria  for  setting  the  coefficient  values  and  a  comparison  of 

the  resulting  performance.  These  include  minimum  mean-square 
f  g  i 

error  3urg  ‘  ,  maximum  likelihood  estimation  Kelley  ancj 

maximum  signal- to-noise  ratio  Shor(4].  it  has  been  shown  that 

all  criteria,  however,  offer  substantial  processino  gain  in 

comparison  with  conventional  delay  and  sum  techniques. 

One  problem  encountered  in  optimum  array  signal  processing 

was  the  requirement  for  knowledge  of  the  noise  correlation 

statistics  in  order  to  calculate  the  filter  coefficient  values. 

This  difficulty  has  been  circumvented,  however,  through  the 

development  of  adaptive  array  processing  methods  by  Widrow[111, 

Chang  and  Tuteur[12],  ComptonU31,  and  Grif  f  i  ths  [  14  1  .  m 

effect,  adaptive  processing  provides  a  method  for  continual 

updating  of  the  filter  coefficients  using  a  simple  feedback  rule. 

Adaptive  processors  may  be  implemented  using  either  analog  or 

diqital  hardware  and  offer  the  additional  advantage  of  being 

able  to  track  slowly  varying  noise  sources. 

An  additional  area  which  has  provided  research  results  in 

array  studies  is  in  the  processing  of  seismic  signals.  Seismic 

processing  techniques  have  been  developed  at  the  MIT  Lincoln 

Laboratory  under  the  Large  Aperture  Seismic  Array  program 

(Capon,  et  al[15])  and  by  several  researchers  concerned  with 

oil  exploration  (Burg^,  Robinson  Claebout  ,  and 

[18] 

Treitel  ) „  Results  in  this  latter  area  are  particularly 
significant  in  that  seismic  "noise"  is  generated  by  the  signal 
source  energy  reflecting  and  reverberating  from  undesired 
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boundaries.  Noise  of  this  type  is  identical  to  the  clutter 

returns  observed  in  the  operation  of  electromagnetic  radar 
systems . 

In  summary,  during  the  past  ten  years  a  large  number  of 
signal  processing  techniques  have  been  developed  for  use  in 
sensor  arrays.  These  methods  provide  increased  spatial  resolu¬ 
tion  and  improved  output  signal- to-noi se  ratio  through  clutter 
and  interference  rejection.  The  purpose  of  the  present  re¬ 
search  program  was  to  demonstrate  that  such  methods  are  directly 
applicable  to  HF  systems  and  to  measure  the  extent  of  their 
improvement. 

D.  APPROACH  TAKEN  IN  PRESENT  STUDY 

A  conventional  receiver  array  beamforming  system  which 
forms  a  single  output  beam  uses  time  delay  steering  as  shown  in 
Fiq.  2.  The  delays  are  arranged  such  that  when  a  planar  wave- 
front  is  incident  on  the  array  from  the  direction  of  interest, 

1<:  tippedrs  ln  Phase  at  the  delay  outputs.  Amplitude  weights 
(Tchel sychev,  raised  cosine,  etc.)  are  added  to  provide  lower 
average  sidelobe  level.  The  beam  output  is  formed  by  summing 
the  delayed  and  weighted  element  signals.  This  output  is 
demodulated  using  a  single-channel  receiver  and  is  then 
processed  to  provide  signal  parameter  information. 

Existing  HF  systems  have  employed  digital  post-receiver 
processing  in  which  an  analog  to  digital  converter  (ADC)  is 
placed  at  the  receiver  output.  Such  a  system  offers  great 
flexibility  in  signal  processing  and  generating  display  formats. 
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The  digital  portions:  of  this  processor  arc  indicated  in  Fia.  2 
usinq  shaded  blocks  and  large  signal  flow  arrows. 

The  array  processino  conf iquration  of  interest  in  this 
research  is  shown  in  Fiq.  3.  Each  antenna  element  is  processed 
usinq  a  tapped-delay-line  filter  of  the  tyj>e  shown  in  Fiq.  1. 

The  coefficients  in  the  delay  lines  are  adjusted  periodically 
usinq  an  alqori  ’  m  which  is  computed  in  the  adaptive  processor. 
The  combination  of  tapped-delay  lines  and  adaptive  algorithm  may 
be  viewed  as  a  time-varying  filter  in  space  and  frequency.  In 
effect,  correct  choice  of  the  coefficients  will  ensure  the 
formation  of  maximum  array  sensitivity  in  the  direction  of 
interest  and  the  formation  of  minimum  sensitivity  --  i.e. 
pattern  nulls  —  in  the  direction  of  interfering  noise  sources. 
The  specification  of  desired  look  direction  and  frequency  struc¬ 
ture  of  the  desired  siqnal  is  made  usinq  reference  inputs  as 
shown  in  Fig.  3.  It  should  be  noted  that  the  connection  between 
the  output  display  and  adaptive  processor  shown  in  this  figure  is 
included  for  purposes  of  generality.  The  alqori thms  reported 
herein  use  only  the  received,  demodulated  signal  and  the  signals 
at  the  adaptive  weights. 

The  digital  processing  blocks  shown  as  shaded  are  used  to 
indicate  one  structure  which  may  be  tsod  to  implement  a  real¬ 
time  adaptive  HF  array  processing  system.  The  delay  lines, 
summinq  operation  and  receiver  in  this  system  are  all  analoq 
devices.  An  alternative,  all  digital  approach  is  shown  in 
Fiq.  4.  In  this  system,  a  phase-coherent  receiver  is  used  to 
demodulate  each  element  output  which  is  in  turn  dioitized  and 
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ref.  inputs 


FIG.  4.  All  digital  adaptive  receiving 
array  beamforming  system 
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processed  using  a  general  purpose  computer.  While  this  system 
is  considerably  more  flexible  than  that  shown  in  Fig.  3,  it 


requires  a  digital  processing  speed  N  times  greater  than  that 
shown  in  Fig.  2,  for  an  array  of  N  elements.  In  choosing 
between  these  alternatives,  parameters  such  as  siqnal  bandwidth, 
number  of  simultaneous  beams,  and  number  of  antenna  elements 


must  all  be  taken  into  account.  This  subject  is  beyond  the 
scone  of  the  present  research  and  the  reader  is  referred  to 
Ref.  19  for  a  more  detailed  discussion.  Althouqh  the  data  used 
m  this  research  program  wore  obtained  by  recording  digitized 
element  output  signals,  the  results  obtained  and  conclusions 
drawn  are  applicable  to  either  the  hybrid  or  all-digital  system 
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III.  EXPERIMENTAL  FACILITIES  AND  DATA  BASE 

A.  EQUIPMENT 

The  data  processed  in  this  report  were  taken  using  CW 
HF  signals  propagated  over  a  2600  Km  east-west  path  between 
Bearden,  Arkansas  and  Los  Banos,  California.  The  equipment  used 
to  collect  these  data  is  part  of  the  wide-aperture  high- 
frequency  radio  research  facility  (WARF)  operated  by  the 
Ionospheric  Dynamics  Group  at  Stanford  Research  Institute,  Menlo 
Park,  California.  A  complete  description  of  this  facility  may 
be  found  in  Reference  20.  The  transmitter,  located  in  Bearden, 
had  a  power  output  of  1  kw  and  was  tuneable  over  the  3-30  MHz  HF 
frequency  band.  The  transmitting  antenna  consisted  of  a  hori¬ 
zontal  log  periodic  antenna  mounted  on  a  100  ft.  tower.  In 
addition,  an  FMCW  transmitter  was  also  available  in  Bearden 
for  purposes  of  obtaining  high  resolution  oblique  ionograms  over 
the  path. 

The  Los  Banos  receiving  array  consisted  of  eight  vertical 
whip  antennas  linearly  spaced  at  a  horizontal  separation  of  360 
meters,  thus  providing  a  total  array  aperture  of  2.5  km.  The 
array  element  output  signals  were  demodulated  using  eiqht 
identical  receivers  driven  by  a  common  local  oscillator  (Fig.  5) . 
A  multiplexed  ADC  was  used  to  sample  the  receiver  outputs 
sequentially  and  store  the  resulting  digital  signals  on  magnetic 
tape.  An  interactive  mini  computer  provided  automatic  data 
collection  as  well  as  operator  documentation  of  the  experimental 
parameters.  Operating  personnel  for  the  tests  were  provided  by 
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FIG.  5.  Los  Banos  receiving  array  system 
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Stanford  Research  Institute. 

B .  DATA  FORMAT 

The  data  processed  in  this  research  were  recorded  using  an 
ADC  sampling  rate  of  500  samples  per  second.  Due  to  the  multi¬ 
plexed  sampling  procedure  employed,  this  provided  an  effective 
per  element  sampling  rate  of  62.5  samples  per  second.  A  total 
of  4096  samples  (512)  samples  per  channel)  were  recorded  per 
data  block  representing  eight  seconds  of  received  data.  Each 
experiment  consisted  of  recording  a  total  of  eleven  such  blocks, 
separated  by  one  minute  intervals.  The  first  data  block  re¬ 
corded  was  a  calibration  block  in  which  the  vertical  whips  were 
disconnected  and  a  fixed  frequency  local  signal  generator  was 
connected  in  phase  to  all  eight  receiver  inputs  as  shown  in 
Fig.  5.  This  calibration  allowed  a  measurement  of  the  amplitude 
and  phase  variations  introduced  by  the  receivers.  During  the 
succeeding  ten  data  blocks,  the  elements  were  reconnected  and 
the  Bearden  transmitter  frequency  was  held  constant.  Experiments 
were  carried  out  under  a  wide  variety  of  ionospheric  conditions 
and  Bearden  transmitter  frequencies.  Oblique  ionograms  were 
recorded  at  regular  intervals  during  the  tests  so  that  propaga¬ 
tion  modes  could  be  identified. 

C.  PRELIMINARY  DATA  PROCESSING  AND  SPURIOUS  SIGNAL  REJECTION 

The  data  used  for  the  analysis  described  in  this  report  were 
originally  collected  by  SRI  personnel  for  an  experiment  to 
determine  the  spatial  characteristics  of  ionospherically  reflected 
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signals  For  this  reason,  the  Bearden  transmitting  signal 

frequency  was  carefully  selected  to  avoid  interference  from  other 
spectrum  users.  Typical  per  element  received  signal-to-noise 
ratios  for  these  tests  were  the  order  of  25  db  to  35  db.  It  is 
well  known  that  under  conditions  of  high  signal-to- 

noise  ratio,  the  optimum  array  processor  consists  of  time  delay 
steering  and  summing  the  element  outputs  —  i.e.  conventional 
array  processing.  Although  this  statement  is  correct  only  if 
the  noise  field  is  isotropic,  the  output  SNR  for  the  conven¬ 
tional  processing  is  for  all  practical  purposes  identical  to  that 
for  optimum  processing  whenever  the  per  element  SNR  exceeds 
about  15  db,  regardless  of  the  spatial  structure  of  the  noise 
field.  For  this  reason,  it  was  decided  that  during  the  present 
investigations,  the  Bearden  signal  would  be  considered  to  be  an 
undesired  interfering  signal  or  jammer,  and  optimum  adaptive 
array  processing  techniques  would  be  used  to  extract  a  weak 
signal,  presumed  to  be  incident  on  the  array  from  some  other 
direction  of  arrival.  The  effectiveness  of  this  processing  was 
then  measured  by  comparing  its  performance  with  that  of  a 
conventional  beams teerinq  system  phased  to  receive  the  some 
weak  signal. 

Preliminary  data  processing  was  employed  in  an  attempt  to 
locate  a  weak  signal  in  the  recorded  data.  The  processing 
method,  which  was  identical  to  that  described  in  Ref.  1,  con¬ 
sisted  of  performing  a  two  dimensional  transform  of  the  original 
data.  Figure  6  rummarizes  this  procedure.  The  original  data 
was  arranged  in  a  matrix  in  which  the  rows  contain  the  time 
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Note  that  the  multiplex 


samples  from  the  eight  receivers, 
system  provides  successive  time  sampies  at  all  eight  elements 
prior  to  returning  to  the  first  element  (i.e.  simultaneous  time 
samples  at  the  elements  are  not  available)  .  A  raised  cosine 
weighting  was  applied  in  both  space  and  time  to  improve  sidelobe 
levels  in  the  final  display.  By  transforming  each  row  and,  in 
turn,  each  column  of  this  matrix,  a  frequency-azimuth  matrix  is 
obtained  in  which  each  row  is  the  received  frequency  spectrum  for 
a  particular  azimuth.  Eight  resolvable  azimuths  may  be  obtained 
using  this  procedure  and  the  azimuthal  spacing  is  given  by 
1/(260x6)  radians  for  a  received  signal  wavelength  of  1  meters. 

The  operation  labelled  skew  compensation  includes  a  linear 

phase  term  to  account  for  the  sequential  element  sampling  and  a 

gain  phase  correction  to  eliminate  local  receiver  differences. 

This  latter  correction  causes  the  calibration  block  data  to 

appear  as  an  ideal  signal  in  the  final  display,  as  shown  in 

Fig.  7.  Similar  displays  for  the  first  and  tenth  data  blocks 

immediately  following  this  calibration  are  shown  in  Figures  6 

and  9.  All  plots  clearly  demonstrate  the  high  signal-to-noise 

ration  of  the  recorded  data.  (Note  that  the  amplitude  scales 
are  linear.) 

In  an  attempt  to  locate  weak  signals  in  these  data,  output 
plots  were  computed  at  gains  of  S,  10,  50  and  100  with  clipping 
applied  to  the  largest  signals.  After  processing  several  blocks 
Of  data  in  this  manner,  a  repeatable  pattern  of  low  amplitude 
energy  was  observed  at  several  points  in  the  two-dimensional 
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Two-dimensional  transform  of  Bearden  signal  data 
for  sixth  data  block  03-04-70  at  0617Z  (5.267  MHz) 


transform  domain.  The  fact  that  this  pattern  was  also  observed 
in  the  calibration  block  indicated  that  the  energy  was  due  to 
spurious  signals  present  in  the  receiving  equipment.  Examples 


of  these  signals  as  found  in  the  data  are  shown  in  Figure  10  and 


11  which  are  Figure  7  and  8  reploted  at  a  gain  of  10. 

An  effective  procedure  was  found  for  eliminating  these  spur¬ 
ious  signals.  First,  the  portion  of  the  transformed  calibration 
block  which  contains  the  true  calibration  signal  energy  is  set 
to  zero.  This  is  a  simple  procedure  because  the  calibration 
frequency  is  known  exactly  and,  as  can  be  seen  from  Figure  7,  it 
occupies  a  relatively  small  portion  of  the  transform  domain. 

The  remaining  energy,  which  consists  of  quantization  noise  and 
spurious  signals,  is  then  subtracted,  point  for  point,  from  the 
transforms  of  the  data  blocks  containing  received  signals,  since 
each  point  of  a  two-dimensional  transform  is  a  complex  number, 
attention  must  be  given  to  both  magnitude  and  phase  in  the  sub¬ 
traction  process.  It  is  reasonable  to  assume  that  the  magnitude 
of  spurious  signals  will  remain  relatively  constant  from  one  data 
block  to  the  next.  The  phase,  however,  will  be  highly  variable 
depending  upon  the  absolute  time  at  which  samples  are  taken. 
Therefore,  to  effectively  eliminate  the  spurious  signal  pattern 
in  the  data,  it  is  sufficient  to  subtract  only  the  amplitude  of 
the  modified  calibration  transform  and  to  leave  the  phase 


transform  unaltered. 


The  effectiveness  of  this  subtraction  procedure  is  shown  by 
comparing  Figure  12,  which  has  been  processed  to  eliminate 
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FIG.  10.  Replot  of  Fig.  7  at  times  ten  gain 


FIG.  11.  Replot  of  Fig.  8  at  times  ten  gain 


spurious  signals  with  Figure  11,  which  has  not.  After  processing 
records  in  this  manner,  gain  factors  of  up  to  100  were  success¬ 
fully  applied  to  locate  weak  signals,  a  quantitative  measure  was 
also  obtained  by  identifying  6  regions  of  particularly  high 
spurious  level  in  the  original  data,  as  shown  in  Figure  13.  The 
reduction  in  spurious  level  for  each  of  these  regions  as  a  func¬ 
tion  of  data  block  number  is  summarized  in  Table  I 


TABLE  1  SUMMARY  OF  SPURIOUS  SIGNAL  REDUCTION  IN  DB 


SPURIOUS 

SIGNALS* 

DATA 
BLOCK  1 

DATA 
BLOCK  2 

DATA 
BLOCK  3 

DATA 
BLOCK  4 

DATA 
BLOCK  5 

1 

-25. Odb 

-19.9 

-21.5 

-24.1 

-16.5 

2 

-30.9 

-11.3 

-15.1 

-16.0 

-16.7 

3 

-19.3 

-12.2 

-15.0 

-14.4 

-17.8 

4 

-  6.0 

-  6.6 

-  5.6 

-11.6 

-  7.9 

5 

-22.2 

-  8.8 

-  9.5 

-10.9 

-16.0 

6 

-14.6 

-  6.1 

-  7.3 

-  7.3 

-13.2 

*As  indicated  on  Figure  13. 


D.  SELECTION  OF  DESIRED  WEAK  SIGNALS 

The  methods  described  in  the  previous  section  were  used 
to  preprocess  several  experimental  data  runs  in  an  attempt  to 
locate  a  weak  signal  in  the  background  noise.  In  each  case, 
spurious-free  transform  plots  similar  to  that  shown  in  Figure  12 
were  manually  scanned.  The  data  selected  for  further  processing 
were  recorded  at  2200  GMT  on  March  3,  1970  at  an  operating 
frequency  of  11.70  MHz.  Figure  14  and  15  show  the  transforms 
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transform  for  block  no.  1  at 
212Z,  03-02-70,  and  frequency  11.7  MHz 
showing  identification  of  desired  signals. 


FIG.  15.  Replot  of  Fiq.  14  at  tines  ten  gain. 
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for  the  first  g  .ta  block  at  unity  and  times  five  gain.  A  total 
of  ten  successive  data  blocks,  taken  at  one  minute  intervals, 
were  available,  in  addition  to  the  initial  calibration  record. 
The  transforms  for  blocks  3,  5,  7  and  9  are  shown  in 
Figures  16  -  19.  All  blocks  indicated  the  presence  of  a  weak 
out  of  band  signal,  denoted  by  d1  in  Figure  14  and  15,  which 

appeared  to  have  fading-like  amplitude  variations  between 
successive  blocks. 

This  signal  was  selected  as  the  first  desired  signal  for 
the  adaptive  array  processor.  It  is  located  at  a  frequency 
offset  of  about  25  Hz  and  is  displaced  in  azimuth  0.51  degrees 
from  the  Bearden  transmitter.  A  second,  in  band,  signal  was 
assumed  to  be  present  in  the  region  denoted  by  d2  in  Figures 
14  and  15.  Due  to  the  presence  of  large  in-band  noise  in  this 
region,  however,  no  direct  evidence  of  the  existence  of  energy 
not  attributable  to  the  Bearden  transmitter  was  ever  obtained. 
The  assumed  azimuthal  deviation  of  this  signal  from  the  nominal 
Bearden  azimuth  was  1.02  degrees. 

The  azimuthal  spread  of  the  Bearden  signal  shown  in 
Figures  14  -  19  is  a  consequence  of  the  large  number  of  iono¬ 
spheric  modes  propagating  during  the  experiment.  An  oblique 
ionogram  which  was  recorded  at  2130Z,  immediately  prior  to  the 
tests,  is  shown  in  Figure  20.  At  least  four  distinct  modes  are 
visible  at  the  transmitter  frequency  (see  arrow) .  Ionograms 
taken  during  the  time  periods  one  hour  preceding  and  following 
-U30Z  indicate  that  the  ionosphere  was  stable  for  the  duration 
of  the  ten  data  blocks  of  interest. 
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It  is  to  be  noted  that  the  rather  extensive  data  processing 
procedure  used  in  this  research  to  identify  the  desired  signals 
dl  and  d2  shown  in  Figure  14  and  15  is  not  normally  included 
as  part  of  an  adaptive  HF  processing  system.  It  was  used  in  the 
present  study  in  order  to  make  most  efficient  use  of  a  data  set 
which  was  recorded  for  other  purposes.  That  is,  the  use  of  two 
dimensional  transforms  f acilitafc-.d  the  identification  of  weak 
signals  located  in  a  strong  interference  background.  In  an 
operational  adaptive  system,  the  direction  and  frequency  struc¬ 
ture  of  these  signals  are  presumed  to  be  known  a  priori  and  this 
information  is  encoded  directly  into  the  adaptive  algorithm,  as 
described  in  the  following  section.  The  noise  field  structure, 
however,  is  assumed  to  be  completely  unknown  a  priori.  The 
motivation  for  the  present  research  is  based  on  the  assumption 
that  real  time  HF  array  processors  do  not  employ  two-dimensional 
processing  techniques  due  to  the  extremely  large  number  of  cal¬ 
culations  required  in  such  a  system. 
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IV.  ADAPTIVE  SIGNAL  PROCESSING  THEORY 
The  theory  of  adaptive  signal  processing  has  been  well 
developed  in  the  last  ten  years  and  applications  have  been 
proposed  for  many  diverse  areas.  The  purpose  of  this  section 
is  to  summarize  the  theory  pertaining  to  the  particular  adap¬ 
tive  algorithm  used  during  the  course  of  the  present  research 
program.  This  algorithm  was  first  proposed  by  Grif f iths . 

A.  MODEL  FORMULATION 

The  antenna  array  of  interest  is  assumed  to  contain  K 
receiving  elements  at  known  spatial  locations.  A  tapped-delay- 
line  processor  (Fig.  1)  is  used  to  process  the  element  outputs. 
Each  line  consists  of  L-l  time  delays  of  A  seconds  each  and  L 
multi  dying  weights. 

The  tapped-delay  line  at  each  antenna  element  is  a  variable 
filter  in  which  the  frequency  shaping  characteristics  of  the 
filter  are  controlled  by  the  tap  settings.  To  illustrate  the 
behavior  of  this  filter,  consider  the  case  of  a  single  delay 
and  two  weights,  as  shown  in  Fig.  21a.  Given  an  input  phasor 
jwot  u 

e  ,  where  ojq  is  a  frequency  such  that  the  delay  is  one 

quarter  wavelength  --  i.e.  A  -  v/2^q,  then  the  output  phasor  is 
j(W0t+TT/2) 

wle  +  w2e  (Fig.  21b).  Assuming  that  the  weights 

can  take  on  both  positive  and  negative  values,  the  gain  and 

phase  response  of  the  filter  w^  +  jW2  can  assume  any  complex 
value . 
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A  multiple  tap  delay  line  (Fig.  22a)  can  be  described  in 
similar  terms.  Weights  spaced  by  A  seconds  are  used  to  set 
the  response  at  frequency  f1  =  1/4A  .  Those  at  a  spacing  of 
2A  seconds  control  the  gain  and  phase  at  one  half  this  frequency, 
f 2  =  1/8A*  Thus,  the  frequency  response  of  a  delay  line  with  L 
weights  can  be  specified  at  L-l  independent,  equally  spaced 
points  between  ^  =  1/4A  and  f^  =  1/4  (L-l)  A  as  illustrated 
in  Fig.  22b.  It  should  be  noted  that  this  discussion  is 
presented  to  illustrate  the  frequency  flexibility  of  a  tapped 
delay  line  processor  and  that  in  an  adaptive  processor  the  tap 
values  are  selected  on  the  basis  of  time  domain  signals  rather 
than  using  a  frequency  spectrum  specification  (see  Section  B 
below).  This  discussion  does,  however,  indicate  the  amount  of 
delay  needed  for  an  array  processor  which  uses  tapped-delay 
line  processors.  Assuming  that  signals  are  received  over  a 
bandwidth  of  B  Hz,  a  minimum  delay  of  A  =  1/2B  seconds  is 
required  to  achieve  unambiguous  frequency  control  over  this 
bandwidth.  If  the  line  contains  L  taps  or,  equivalently,  a 
total  delay  of  (L-1)/2B  seconds,  the  minimum  frequency  resolution 
will  be  B/L  Hz.  For  example,  an  HF  system  operating  at  a 
100  kHz  bandwidth  which  uses  a  five  tap  delay  line  would  be 
capable  of  frequency  discrimination  within  20  kHz  bandwidths 
and  would  require  individual  delays  of  5  ysec.  The  total  delay 
line  length  would  be  20psec. 

The  array  output  signal  for  the  processor  shown  in  Fig.  1 
is  formed  from  the  summation  of  all  delayed  and  weighted  signals. 
Denoting  the  set  of  KL  weights  in  this  processor  by  a  vector  W 
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and  the  set  of  KL  signals  present  at  each  weight  at  time  t  by 
the  vector  X(t) 


X(t) 


/ 


(t) 

x2(t) 


(2) 


XKL^ 

then  the  array  output  signal  y(t)  can  be  expressed  as  an 
inner  product,  i.e. 


y(t)  =  WTX(t) 


KL 


l  w-x.ft). 


i=l 


(3) 


This  notation  may  be  easily  modified  to  include  the  effects  of 
adaptation  in  which  the  value  of  each  weight  w±  changes  at 
each  adaptive  step.  Let  W(k)  express  the  value  of  the  weights 
after  k  adaptations,  then 


W  (k)  = 


j  wx(k) 

w2  (k) 


\wKL(k) 


t-hs  output  signal  at  time  t  becomes 


T 


y(t)  =  w1  (k)x(t) 


(4) 


(5) 
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It  is  presumed  that  the  k^1  adaptation  occurs  at  time  t.  . 

K 

The  array  processor  studied  in  this  research  is  that  which 
minimizes  the  mean-square  error  between  a  desired  signal  wave¬ 
form  d(t)  and  the  array  output  signal  y(t).  The  desired 
waveform  was  assumed  to  be  planar  and  to  propagate  across  the 
array  from  a  known  direction  of  arrival  0^.  Each  antenna 
element  is  also  assumed  to  receive  noise  signals  consisting  of 
both  directional  interference  and  non-directional  background 
sky  noise.  Although  the  precise  nature  of  these  noises  — 
i.e.  direction  of  arrival  and  spectrum  —  is  unknown,  they  are 
all  assumed  to  be  uncorrelated  with  the  desired  signal. 

The  set  of  coefficients  W*  which  minimizes  the  mean- 
square  error  between  y(t)  and  d(t)  is  the  Wiener  filter 
which  has  been  discussed  by  many  authora^6' 8,9 ' It  is  given 
by  the  solution  to  the  normal  equations 

-*  =  “XX-X  (6) 

a 

where  gxx  is  the  autocorrelation  matrix  of  the  signals  ob¬ 
served  at  the  weights  and  PY  is  the  cross-correlation  vector 

—xd 

between  these  signals  and  the  desired  signal.  That  is 

Rxx  =  E[X(t)XT(t)]  (7) 

Px^  =  E  [d (t)  X(t)  ]  (8) 

where  E [ - ]  denotes  expectation,  or  time  average,  and  superscript 
T  denotes  transpose. 

In  general,  an  array  processor  similar  to  that  shown  in 
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Fig.  1  with  the  tapped-delay  line  weights  set  to  the  values 

given  in  Eq.  (6)  will  have  a  mainlobe  in  the  direction  6  of 

d 

the  desired  signal  and  will  also  contain  pattern  nulls  in  the 
direction  of  strong  interfering  noise  sources.  Clearly,  a 
processor  of  this  type  can  achieve  considerable  improvement  in 
output  signal-to-noise  ratio  as  compared  with  conventional 
time-delay  and  sum  processing  under  conditions  of  low  per  element 
signal-to-noise  ratio.  Optimization  criteria  other  than  minimum 
mean-square  error  have  been  studied  in  the  literature— e .g.  see 
References  4,  13  and  15  —  and  these  processors  also  achieve 
desirable  antenna  array  patterns.  The  reason  for  choosing  a 
mimmum-mean-square  error  criterion  in  the  present  research  was 
that  a  simple  adaptive  algorithm  is  available  which  achieves 
this  solution  using  a  minimal  number  of  computations.  Similar 
algorithms  for  other  performance  criteria  have  not  yet  been 
developed. 

B  THE  ADAPTIVE  ALGORITHM 

Due  to  a  lack  of  knowledge  regarding  the  noise  structure, 
the  optimum  weight  set  W*  defined  by  Eq.  (6)  cannot  be  cal¬ 
culated  directly.  Procedures  are  available tl5J  for  measuring 
the  correlation  statistics  at  the  element  outputs  and  using  this 
information  to  compute  the  matrix  and  the  resulting  optimum 

weight  set.  One  difficulty  with  this  method,  however,  is  that 
the  measurement  of  correlation  parameters  and  weight  set  compu¬ 
tation  is  complex,  requiring  numerous  numerical  calculations. 

In  addition,  the  weight  set  so  obtained  will  be  incorrect  in 


the  presence  of  chanqinq  noise  fields. 

An  alternative  procedure  for  computing  W*  is  the  adaptive 

r  i4i 

algorithm  first  proposed  by  Griffiths1  .  In  this  method,  the 
weight  vector  is  recalculated  at  successive  times 
t^, t2, . . .t^, . . .  using  the  equation 

w(k)  =  w(k-l)  +  u c£x  -  y(tk)x(tk)],  (9) 

d 


where:  tk  =  time  at  which  k*"*1  adaptation  takes  place 

W(k-l)  =  set  of  filter  coefficients  used  prior  to  tk 
W(k)  *  set  of  filter  coefficients  used  after  tk 
y(tk)  =  array  output  signal  at  time  tk 
X(tk)  =  vector  of  signals  stored  in  the  filter 
V  =  proportionality  constant  which  controls  the 
rate  of  convergence  of  the  algorithm 

A 

P  =  estimate  of  signal  cross-correlation  vector 
~*d 

A  complete  description  of  the  derivation  and  convergence  pro¬ 
perties  of  this  algorithm  may  be  found  in  Reference  14.  It  can 
be  shown  that  for  a  value  of  p  which  satisfies 

KL  , 

0  <  p  <  2/  l  E{xf(t.)],  (10) 

i=l  1  * 

KL  2 

where  l  E(xf(t.  )]  is  the  total  signal  power  present  in  the 
i=l  1  * 

array  processor,  then  W(k)  converges  in  the  mean  to  W*  as  k 
tends  to  infinity.  That  is,  the  limiting  mean  value  of  W(k)  is 
equal  to  W*  and  the  standard  deviation  is  a  non-zero  value 
which  is  proportional  to  the  constant  p. 


Several  practical  factors  must  be  considered  when  Eq.  (9) 
is  implemented  as  an  array  processor.  The  choice  of  the  optimum 
value  for  y  involves  a  tradeoff  between  the  convergence  rate 
and  the  steady-state  adaptation  output  noise  which  is  produced 
by  the  fact  that  the  weights  are  constantly  changing.  A  value 
near  the  upper  limit  in  Eq.  (10)  will  produce  a  rapid  conver¬ 
gence  rate  and  relatively  large  output  adaptive  noise.  Values 
near  zero,  however,  have  the  opposite  effect  —  i.e.  long  con¬ 
vergence  times  and  negligible  output  noise.  Given  that  y  must 
be  large  enough  to  permit  the  processor  to  track  time  variations 
in  the  received  noise  field,  an  effective  value  can  be  deter¬ 
mined.  As  shown  in  Section  V,  the  convergence  properties  are 
relatively  insensitive  to  small  changes  in  y  and  values  between 

.04  and  .40  of  the  upper  limit  in  (10)  produced  satisfactory 
results . 

An  additional  factor  which  must  be  considered  is  the 

estimated  cross-correlation  vector  P  .  Given  that  the  desired 

d 

signal  d(t)  is  assumed  to  be  planar  and  incident  on  the  array 
from  a  known  direction  6d,  the  value  of  this  vector  may  be 
calculated  from  a  knowledge  of  the  power  spectrum  of  d(t) .  For 

the  linear  array  geometry  of  interest  in  this  research,  P 

.  ~Xd 

has  a  particularly  simple  form,  if  At  is  the  interelement 

delay  for  signals  arriving  from  direction  ed,  the  ith 

component  of  this  vector  (which  is  the  cross-correlation  between 

the  signal  observed  at  the  ith  filter  coefficient  and  the  desired 
array  output  signal)  is  given  by 
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(11) 


Pi  =  rd(To  -  nAr  -  mAT)  . 

In  this  equation,  rd(*)  autocorrelation  function  of  the 

desired  signal  which  is  computed  as  the  inverse  of  the  given 

power  spectrum,  n  is  the  receiving  element  number  and  m  is 

the  number  of  delays  preceding  the  coefficient  of  interest.  The 

symbol  tq  is  the  reference  delay  of  the  desired  signal  which 

is  generally  set  equal  to  the  delay  to  the  center  tap  in  the 

delay  line  connected  to  a  central  element  in  the  array.  Thus, 

the  desired  array  output  signal  is  assumed  to  be  d(t-x  )  when 

o 

d ( t)  is  the  signal  received  by  the  element  nearest  the  desired 
source.  The  performance  of  the  array  processor  can  be  shown 
to  be  relatively  insensitive  to  changes  in  tq  so  long  as 

a  value  near  that  corresponding  to  the  mid-point  delay  of  the 
processor  is  used. 
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v.  RESULTS 

This  section  summarises  results  obtained  from  the  off-line 
processing  of  ten  successive  data  blocks  taken  over  the  Bearden- 
Los  Banos  path,  as  described  in  Section  III.  The  data  were 
processed  to  demonstrate  the  advantages  of  optimum  array  pro¬ 
cessing  techniques  in  low  signal-to-noise  ratio  environments. 

TWO  desired  signals  of  known  center  frequency  and  direction  of 
arrival  were  assumed  to  be  present  in  the  data  -  see  Pigs.  14 

and  15 '  ^  BSarden  Si9na1'  «hioh  constituted  the  major  faction 

of  the  received  energy,  was  assumed  to  be  an  undesired  inter¬ 
fering  signal  of  unknown  spectral  and  spatial  content.  Least 
mean-square  error  array  processors  were  developed  to  extract 
each  of  the  desired  signals,  in  turn.  The  performance  of  these 
processors  was  then  compared  with  conventional  time  shift  and 

spatial  weighting  methods  of  receiveing  the  same  desired  signals. 
A.  OPTIMUM  PROCESSOR  CALCULATIONS 

As  described  in  Section  IV  above,  the  set  of  filter 
coefficients  which  minimises  the  mean-square-error  perfor¬ 
mance  criterion  can  be  calculated  from  a  knowledge  of  the  auto¬ 
correlation  matrix  of  the  received  signals  and  the  estimated 

cross-correlation  vector  for  the  desire 

me  desired  signal  (see  Eqs .  (6), 

(7) , (8) ) .  The  array  performance  achieved  with  this  set  of 
coefficients  is  equal  to  that  achieved  on  the  average  by  an 
adaptive  processor  and  therefore  represents  a  performance  base¬ 
line  which  may  be  used  to  evaluate  the  effectiveness  of  minimum 
mean-square-error  array  processing. 
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The  optimum  filter  coefficients  w*  were  calculated  for 
each  of  the  ten  successive  data  blocks.  Correlation  statistics 
were  computed  by  using  time  averages  over  the  eight  seconds  of 
data  (512  samples  per  receiving  element)  within  each  data  block. 
Statistics  of  this  type  were  measured  for  processors  consisting 
of  two,  five,  and  ten  filter  coefficients  per  array  element. 
Since  these  statistics  depend  only  upon  the  properties  of  the 
received  signals,  they  are  not  dependent  upon  the  direction  of 
arrival  and  frequency  of  the  presumed  desired  signal.  This 
latter  information  is  contained  in  the  structure  of  the 
estimated  cross-correlation  vector  P  defined  in  Eq.  (8). 

a 

The  desired  signals  used  throughout  this  research  were 
presumed  to  have  a  sinusoidal  CW  time  structure  and  planar 
spatial  wavefronts.  It  was  assumed  that  the  phase  of  these 
Signals,  as  measured  at  the  array  element  nearest  the  desired 
source,  could  be  modelled  as  a  random  variable  uniformly  dis¬ 
tributed  over  0  to  2ir  radians.  This  assumption  is  consistent 
with  experimental  observations  of  ionospherically  reflected, 
single  mode  CW  signals.  The  autocorrelation  function  for  CW 
signals  with  random  phase  is  given  by 

a2 

rd(T)  =  T  cos  woT  (12) 

where  A  is  the  received  signal  amplitude  and  u>  is  the 
center  frequency.  The  components  of  the  cross-correlation 

A 

vector  P_x  for  an  array  processor  with  KL  coefficients  are  then 
given  by 


48 


COS  (11  (t  -T.) 
O  O  1 


(13) 


In  this  expression,  is  the  time  delay  between  the  i^ 

filter  coefficient  and  the  element  nearest  the  desired  source 


and  t  is  the  reference  delay  of  the  desired  signal  (see 
Eq.  (11)).  The  method  used  to  calculate  this  delay  was 


where  t  ,  and  t  are  the  minimum  and  maximum  values  of  t- 
min  max  i 

in  Eq.  (13)  . 

The  only  parameter  in  Eq.  (13)  which  cannot  be  calculated 
from  a  knowledge  of  the  direction  of  arrival  and  frequency  of 
the  desired  signal  is  the  amplitude  A.  It  is  not  reasonable 
to  assume  that  this  parameter  is  known  a  priori  in  most  appli¬ 
cations  since  knowledge  of  A  implies  that  the  received  per 
element  signal- to-noise  ratio  is  also  known.  Inspection  of 
Eq.  (6),  however,  shows  that  only  the  gain  of  the  optimum  filter 
W*  changes  with  varying  A.  Because  the  processor  output  signal 
to-noise  ratio  is  independent  of  gain  —  i.e.  increased  gain 
increases  both  the  output  signal  power  and  the  output  noise 
power  —  an  arbitrary  gain  A  *  /2  can  be  used  without  affecting 
the  detection  performance  of  the  processor. 
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B.  CONVENTIONAL  PROCESSOR  CALCULATIONS 

A  second  type  of  fixed,  non-adaptive  processor  was  computed 
and  tested  usinq  the  received  data.  This  processor  consisted  of 
time  delays  and  amplitude  weighting  as  shown  in  Fig.  23.  The 
time  delays  At , 2At , . . . 7At  are  chosen  such  that  a  planar  signal 
arriving  from  the  known  direction  0^  appears  in  phase  at  the 
output  of  the  delays .  Delay  values  were  calculated  for  each  of 
the  two  desired  signals  d^  and  d2  shown  in  Fig.  15.  The 
amplitude  weights  ci'c2'‘"c8  were  selected  using  a  raised 
cosine  in  which  a  ten  element  array  was  assumed  with  the  outer 
two  elements  receiving  zero  weight  values.  The  weights  for  the 
remaining  eight  elements  are  given  in  Table  2. 

TABLE  2 

AMPLITUDE  WEIGHTS  FOR  CONVENTIONAL  PROCESSOR 


Element  No.  Weight 

1  .17 

2  .41 

3  .75 

4  .97 

5  .97 

6  .75 

7  .41 

8  .17 


The  antenna  array  pattern  for  the  conventional  processor 
steered  toward  signal  d^  is  shown  in  Fig.  24.  A  linear 
amplitude  scale  has  been  used  in  this  display  and  the  separation 
between  adjacent  azimuthal  tick  marks  is  0.510  degrees.  Note 
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signal 


FIG.  23.  Non-adaptive  time-delay  and  amplitude 
weighting  beamforming  system 
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that  only  the  mainlobe  pattern  has  been  plotted  to  simplify  the 
display.  Because  of  the  wide  spacing  (360m)  between  the  individ¬ 
ual  antenna  elements,  the  complete  pattern  contains  grating  lobes 
at  increments  of  4.08  degrees  for  signals  received  at  11.70  MHz. 

The  mainlobe  pattern  for  the  conventional  processor  steered  to¬ 
ward  the  second  desired  signal  d2  is  identical  to  that  shown 
in  Fig.  24,  except  that  it  is  shifted  one  increment  to  the  right. 

Antenna  patterns  drawn  for  the  conventional  processors  in  this 
research  are  frequency  independent  over  the  bandwidth  of  the 
received  signals.  This  is  due  to  the  fact  that  the  processors 
are  assumed  to  consist  of  ideal  time  delays  and  multiplying 
weights.  In  practice,  narrow  band  conventional  processors  are 
constructed  using  90°  phase  shift  elements  at  the  array  center 
frequency,  which  will  generally  produce  a  frequency  dependent 
array  pattern.  Due  to  the  relatively  low  sampling  rate  of 
62.5  Hz  used  in  the  present  experiment,  however,  the  received 
signals  were  bandlimited  to  31.25  Hz  and  90°  phase  shifters 
would  effectively  provide  ideal  time  delay  over  this  bandwidth. 

In  summary,  the  conventional  pattern  for  the  signals  recorded 
in  this  work  is  frequency  independent  for  both  time  delay  and 
phase  shift  steering. 

An  alternative  method  for  displaying  antenna  array  patterns 
is  illustrated  in  Figs.  25  and  26  which  show  the  conventional 
patterns  for  d^^  and  d2,  respectively.  This  procedure  pre¬ 
sents  the  array  response  in  the  eight  resolvable  azimuthal 
directions  as  a  function  of  frequency  over  the  31.25  Hz  bandwidth. 
The  ninth  azimuthal  line  shown  in  these  figures  is  a  duplicate 


of  the  first  line  and  is  included  for  purposes  of  symmetry. 

Note  that  these  figures  demonstrate  the  frequency  independence 
of  the  conventional  processor. 

C.  PERFORMANCE  COMPARISONS  FOR  DESIRED  SIGNAL  D^. 

The  previous  discussions  have  described  two  methods  of 
computing  fixed  processors.  The  performance  of  these  processors 
can  be  effectively  compared  by  computing  the  resulting  output 
signal-to-noise  ratio  for  each  processor.  This  calculation  was 
carried  out  by  first  computing  the  total  output  signal  power  in 
terms  of  the  known  signal  parameters  and  the  processor  coeffi¬ 
cients.  The  sum  of  signal  plus  noise  cutput  power  for  the 
recorded  data  was  then  measured  by  implementing  the  processor 
and  time  averaging  the  square  of  the  output  signal.  An  equi¬ 
valent,  direct  calculation  can  be  made  using  the  measured  auto¬ 
correlation  matrix. 

A  summary  of  the  output  signal-to-noise  ratio  performance 
for  both  optimal  and  conventional  processors  designed  to 
receive  signal  dx  is  given  in  Table  3,  as  a  function  of  the 
data  block  number.  The  individual  element  signal-to-noise 
ratios  are  also  included  for  purposes  of  comparison.  Optimal 
filters  consisting  of  two  weights  per  element  (16  weights  total), 
five  weights  per  element  (40  weights  total),  and  ten  weights 
per  element  (80  weights  total)  were  computed  for  the  first  and 
third  data  blocks  and  the  40  weight  processor  was  computed  for 
all  ten  cases.  The  average  array  gain  —  computed  as  the 
improvement  above  single  element  SNR  —  for  these  processors 
is  given  in  Table  4  below. 
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TABLE  4 


AVERAGE  PROCESSOR  GAIN 

PROCESSOR  TYPE  GAIN  OVER  SINGLE  ELEMENT 

conventional  12.5  db 

16  weight  optimal  17.4 

40  weight  optimal  32.5 

80  weight  optimal  41.6 

Figures  28  through  36  illustrate  the  antenna  array  patterns 
achieved  for  the  optimal  processors.  A  comparison  of  the  16, 

40,  and  80  tap  filters  for  block  1  is  shown  in  Figs.  27-29  and 
for  block  3  in  Figs.  30-32.  The  performance  of  the  40  weight 
optimal  filter  on  data  blocks  5,  7,  and  9  is  given  in  Figs.  33-35. 
Finally,  a  conventionally  displayed  array  pattern  is  shown  in 
Fig.  36  for  the  40  weight  processor  computed  for  block  1.  The 
pattern  is  plotted  at  three  frequencies:  the  Bearden  trans¬ 
mitter  frequency  (minimum  response) ,  the  d1  desired  frequency 
(maximum  response)  and  an  intermediate  frequency.  Note  that 
all  of  the  optimal  patterns  have  a  maximum  in  the  direction  of 
the  desired  signal,  but  that  this  maximum  is  frequency  dependent 
and  appears  only  at  the  frequency  of  the  desired  signal.  Most 
patterns  also  show  a  spatially  broad,  deep  null  at  the  Bearden 
frequency. 

It  should  be  noted  that  block-to-block  variations  in  the 
output  SNR  of  the  optimal  filters  shown  in  Table  3  is  due  to  two 
factors:  1)  The  optimal  filter  calculation  depends  on  the 

correlation  statistics  of  the  received  data  which  may  vary  from 
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azimuth  (0.51°  per  increment) 


Single  frequency  antenna  array  patterns  for 
optimal  filter  shown  in  Fig.  28. 


block  to  block  and  2)  The  received  signal-to-noise  ratio 
varies  between  blocks.  An  example  of  the  latter  factor  can  be 
found  by  comparing  the  40  weiqht  optimal  filters  for  blocks 
1  and  2.  The  second  block  has  a  processed  output  SNR  which  is 
13.3  db  higher  than  that  for  the  first  data  block.  Inspection 
of  the  elemental  SNR's,  however,  shows  that  11.4  db  of  this 
change  can  be  attributed  to  a  drop  in  the  received  SNR.  Varia¬ 
tions  due  to  correlation  changes  are  important  when  considering 
the  potential  use  of  an  adaptive  processor  because  they  require 
corresponding  changes  in  the  filter  coefficients.  Received  SNR 
variations  are  relatively  unimportant  to  an  adaptive  processor 
because  they  reflect  simple  gain  changes  and  do  not  require 
adaptation . 

Calculations  were  made  to  determine  the  time  rate  of  change 
of  the  correlation  statistics.  The  procedure  used  was  to  measure 
the  performance  of  an  optimal  filter  computed  for  one  particular 
data  block  by  computing  output  SNR  for  other  data  blocks.  In 
effect,  this  procedure  is  equivalent  to  computing  an  optimal 
filter  during  one  data  block  and  then  fixing  the  weights  for  use 
in  other  data  blocks.  Figure  37  dhows  the  performance  change 
observed  when  the  optimal  40  weight  filter  for  the  first  data 
block  was  applied  to  all  succeeding  blocks.  Similar  plots  for 
optimal  filters  computed  for  blocks  3  and  7  and  applied  to  other 
data  sets  are  shown  in  Fig.  38.  These  plots  indicate  that  the 
penalty  for  not  adapting  to  track  the  optimal  filter  is  the 
order  of  7  to  10  db,  depending  upon  the  length  of  time  between 
the  calculation  of  the  optimal  filter  and  the  SNR  measurement. 
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D.  PERFORMANCE  COMPARISONS  FOR  DESIRED  SIGNAL  D 

2 

Inspection  of  the  two  dimensional  array  patterns  for  the 
optimal  processors  given  in  Figs.  27-36  shows  that  the  major 
improvement  achieved  in  output  SNR  is  due  to  the  null  which  is 
Placed  on  the  Bearden  transmitter.  In  order  to  determine  the 
performance  improvement  of  optimum  processing  for  desired  signals 
which  are  less  isolated  in  the  data  displays  shown  in  Figs.  14-19, 
n  second  desired  signal  d2  was  assumed  to  be  present  at  the 
Bearden  frequency  and  displaced  1.02  degrees  from  the  Bearden 
azimuth  as  shown  in  Fig.  15.  A  time  delay  and  amplitude  shaded 
conventional  processor  was  computed  for  this  signal.  The  re¬ 
sulting  antenna  pattern  is  shown  in  Fig.  26.  The  pattern  for  a 
40  weight  optimal  processor  computed  using  signals  in  the  first 
data  block  is  shown  in  Fig.  39.  Note  that  the  response  has  a 
spatial  null  at  the  Bearden  frequency  for  all  azimuths  except 
that  corresponding  to  the  assumed  direction  of  arrival  of  a,. 

The  per  element  SNR  for  this  case  was  -38.1  db,  conventional2 
processing  produced  an  output  SNR  of  -30.9  db  and  optimal  pro- 
cessinq  resulted  in  a  -25.0  db  output  SNR. 

The  rather  low  resulting  antenna  gain  for  the  optimal 

processor  is  a  consequence  of  the  fact  that  the  Bearden  signal 

is  widely  spread  in  azimuth.  Thus,  at  the  Bearden  frequency, 

the  int, rfering  noise  is  nearly  isotropic  and  conventional 

processing  is  near  optimal.  The  five  db  improvement  achieved  by 

the  optimal  processor  is  due  to  spatial  irregularities  in  the 
Bearden  signal. 
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frequency 


FIG.  39. 


Optimal  filter  with  40  weights  designed  for  data 
block  1  and  desired  signal  d2» 
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E.  ADAPTIVE  ALGORITHM  IMPLEMENTATION 

The  adaptive  algorithm  given  in  Eq.  (9)  was  implemented  on 
the  recorded  data  using  the  out  of  band  siqnal  d^  as  the 
desired  signal.  The  processor  used  in  these  studies  contained 
a  total  of  40  adapted  filter  coefficients  --  i.e.  5  taps  per 
element.  Several  initial  weight  values  and  adaptation  time 
constants  were  studied.  In  all  cases,  the  performance  of  the 
adaptive  processor  showed  excellent  agreement  with  theoretical 
calculations  based  on  the  measured  correlation  statistics  of  the 
data. 

Figure  40  illustrates  the  behavior  of  adaptive  filter  in 
the  region  of  the  optimal  solution.  The  vertical  axis  shows 
the  processor  output  SNR  for  the  optimal  set  of  filter  coeffi¬ 
cients  (-19.42  db)  and  for  4  adaptive  processors  implemented 
using  proportionality  constants  y  =  200,  400,  1000,  and  2000 
[see  Eq.  (9)].  The  data  used  for  these  tests  was  that  recorded 
in  the  first  data  block.  In  all  four  adaptation  experiments, 
the  filter  coefficients  were  initially  set  to  the  optimal  values 
previously  computed  for  this  data  block.  Thus,  output  SNR  of 
the  four  processors  is  equal  to  the  optimum  value  before  adap¬ 
tation  begins.  As  adaptation  progresses,  however,  the  adaptive 
filters  depart  from  the  optimal  \alue  by  an  amount  proportional 
to  y,  the  adaptive  step  size. 

The  suboptimality  of  the  adaptive  processors  shown  in 
Fig.  40  is  a  well  known  consequence  of  the  fact  that  the  weight 
computation  algorithm  given  in  Eq.  (9)  continually  adjusts  the 


70 


filter  coefficients.  Widrow[22]  has  termed  this  behavior 
mis  adjustment  and  has  given  formulas  for  calculating  the  average 
percentage  offset  as  a  function  of  the  proportionality  con¬ 
stant  y.  it  should  be  noted,  however,  that  even  the  largest 
u  value  used  in  these  experiments  produced  a  steady-state 
offset  of  less  than  1  db . 

The  ranee  of  proportionality  constants  used  in  these  experi¬ 
ments,  200  to  2000,  can  be  related  to  the  upper  bound  given  by 
Eq-  (10)  which  is  repeated  below  for  reference, 

KL  o 

0  <  ^  <  2/  I  E[x2(tk)J  .  (10) 

The  bound  can  also  be  computed  as  two  over  the  sum  of  the 

diagonal  elements  in  the  matrix.  For  the  first  data  block, 

the  value  is 

0  <  P  <  5,300  .  (15) 

Thus,  the  values  ranged  between  0.038  and  0.38  of  the  upper 
1  im  i  t . 

A  series  of  adaptive  experiments  were  also  carried  out  to 
determine  the  time  constant  of  adaptation,  as  a  function  of  the 
proportionality  constants.  in  these  experiments,  the  initial 
values  of  the  filter  coefficients  were  chosen  to  be  far  removed 
from  the  optimal  values  and  the  relaxation  time  was  measured 
as  the  filter  adapted  toward  the  optimum  solution. 

Figures  41-44  illustrate  the  results  of  a  time  constant 
determination  experiment  using  data  block  number  4.  In  each 
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fig.  42.  Adapted  output  SNR  with  trend  1 
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No.  of  adaptations 


FIG.  4 


Adapted  output  SNR  with  trend  1 
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w-ts  initially  not  to  a  random  value  boforo 
adaptation  began.  These  figures  show  the  output  SNR  as  a  function 
of  number  of  adaptations  for  values  of  p  of  200,  400,  1000  and 
2000.  A  best-fit  linear  trend  line  was  also  computed  and  plotted 
for  each  experiment.  Inspection  of  these  results  shows  that 
adaptation  time  constant  decreases  with  increasing  p,  resulting 
in  a  faster  adaptation  rate,  and  that  the  adaptive  noise  --  i.e. 
variations  about  the  trend  line  -  increases  with  increasing  p. 
Both  of  these  observations  agree  with  theoretical  predictions 
(see  Ref  s  [14]  and  [22]).  Table  5  below  summarizes  the  time 
constant  measurements  taken  from  the  trend  line  on  these  curves. 


a  nu  j-jaj  j 


ADAPTIVE  TIME  CONSTANT  MEASUREMENTS 


PROPORTIONALITY 
CONSTANT  p 


ADAPTATION 

RATE 


NO.  OF  ADAPTATIONS 
TO  HALF  POWER 


200 

0.85*10-3db/adapt 

3,530 

400 

1. llxio-3 

2,700 

1000 

1. 73><10"3 

1,730 

2000 

4.38xl0-3 

680 

Complete  theoretical  predictions  of  adaptive  time  constants 
have  been  carried  out123'24’  in  terns  of  the  correlation  statis- 
tics  of  the  data.  It  can  be  shown  that  the  adapted  mean  square 
error  converges  as  a  sum  of  exponentials  in  which  the  exponent 
rate  depends  on  the  eigenvalues  of  the  £xx  correlation  matrix 
and  on  the  value  of  the  proportionality  constant.  Specifically, 
it  can  be  shown t24]  that  the  adapted  mean  square  error  after  k 
adaptations  e^(k)  is  given  by 
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steady-state  adapted  mean-square  error 

constant  depending  upon  initial  set  of 
coefficients  (prior  to  adaptation) 

(1  -  yXi) 

. th  .  ,  _ 

l  eigenvalue  of  R  . 

:AA 

Because  of  the  constraints  on  \i,  each  ri  has  a  magnitude 
value  less  than  unity.  Thus,  the  second  term  in  (16),  which 
is  the  transient  term,  tends  to  zero  as  k  -*  «>.  It  should  be 
noted  that  the  value  of  e^s  depends  upon  y  and  is  somewhat 
hiqher  than  the  optimum  value  as  shown  experimentally  in 
Fig.  41. 

Although  the  convergence  rate  of  a  sum  of  exponentials 
cannot  be  accurately  approximated  by  one  exponential  rate,  it 
can  be  bounded  in  terms  of  the  largest  and  smallest  exponent. 
The  fastest  and  slowest  possible  convergence  rates  are  given  by 


where , 
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where  rmax  and  rmin  are  the  maximum  and  minimum  values  of 
r i >  respectively.  These  values  were  computed  for  the  fourth 
data  block  in  terms  of  the  eigenvalues  of  the  R  matrix  given 

A  A 
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in  Table  6. 


TABLE  6 

EIGENVALUES  OF  Rxx  MATRIX  FOR  4th  DATA  BLOCK 


NO. 

EIGENVALUE 

1 

7.22798E-05 

2 

3. 13664E-05 

3 

6.6189  8E-06 

4 

4. 41183E-06 

5 

1. 31422E-0  6 

6 

5 . 730  37E-0  7 

7 

4. 72674E-07 

8 

3.69323E-07 

9 

2. 57834E-07 

10 

1. 89  351E-0  7 

11 

1. 54275E-07 

12 

1. 07782E-07 

13 

7 . 24330E-0  8 

14 

4 . 4210  7E-0  8 

15 

3. 68840E-08 

16 

2 . 839 14E-0  8 

17 

2 . 11313E-0  8 

18 

1. 65635E-0  8 

19 

1 . 21556E-0  8 

20 

7.77124E-09 

NO. 

EIGENVALUE 

21 

6 . 77182E-05 

22 

2.98098E-05 

23 

6 . 05305E-06 

24 

4.216  84E-06 

25 

1 . 04394E-06 

26 

5 . 03522E-0  7 

27 

3.99563E-07 

28 

3. 39554E-07 

29 

2 . 17305E-07 

30 

1 . 6  32  32E-0  7 

31 

1. 30565E-07 

32 

8 . 5060 1E-0  8 

33 

6 . 19560E-0  8 

34 

4 . 02069E-0  8 

35 

3. 05655E-0  8 

36 

2 . 7  3219E-0  8 

37 

1.82156E-08 

38 

1. 56636E-08 

39 

8. 17442E-09 

40 

5 . 54  834E-09 

Assuming  that  the  output  SNR  convergence  rate  is  similar  to  that 
given  for  the  mean-square  error  in  (16),  the  adaptation  rates 
can  be  bounded  using  the  maximum  and  minimum  eigenvalues  in 
Table  6.  That  is, 


r  =  1  - 

max 

yx7. 23><10"5 

(19) 

r  .  =  1  - 

min 

Ux5.55xlo-5 

(20) 

Table  7  shows  the  results  of  these  calculations  which  indicate 

that  the  observed  time  constants  are  well  within  the  theoretical 
bounds  . 

78 


■'  .  . . - 


TABLE  7 


THEORETICAL  BOUNDS  ON  ADAPTATION  RATES  FOR  4th  DATA  BLOCK 


PROPORTIONALITY 

MAXIMUM 

MINIMUM 

CONSTANT  y 

THEO.  RATE 

THEO.  RATE 

200 

0.126  db/adapt 

0.97xlo“5  db/adapt 

400 

0.161 

1.93*10"5 

1000 

0.653 

4.78xl0"5 

2000 

1.360 

9.65xl0~5 
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VI.  CONCLUSIONS 


The  research  results  described  in  this  report  were  obtained 
by  processing  digital  recordings  of  wide  aperture  HF  array  data 
taken  at  the  WARF  facility  in  Los  Banos,  California.  These 
recordings  were  taken  early  in  1970  over  a  one-hop,  forward 
scatter  ionospheric  path  using  a  transmitter  located  in  Bearden, 
Arkansas.  The  data  were  originally  collected  for  an  experiment 
concerned  with  measuring  the  ultimate  spatial  resolution  of  a 
large  HF  array.  In  effect,  the  prrpose  cf  the  experiment  was 
to  determine  the  degree  to  which  ionospheric  irregularities 
limit  the  spatial  coherence  of  one-hop  reflected  signals.  For 
this  reason  the  signals  were  recorded  under  extremely  high 
signal- to-noise  ratio  conditions.  The  operating  frequency  was 
carefully  selected  to  avoid  interference  and,  because  iono¬ 
spheric  doppler  shifts  are  seldom  more  than  a  few  Hertz,  a  very 
narrow  receiver  bandwidth  (31.25  Hz)  was  used. 

Receiving  parameters  such  as  these  are  not  typical  of 
most  HF  systems  which  utilize  bandwidths  as  large  as  100  KHz 
and  are  often  interference  limited.  However,  the  data  were 
recorded  at  eight  spatially  separated  elements  and  thus  could 
be  used  to  examine  the  potential  advantages  of  adaptive  beam- 
forminq  procedures.  The  conclusions  drawn  from  these  measure¬ 
ments  were  encouraging  but  must  necessarily  be  restricted  to 
narrow-band  systems  operating  in  the  presence  of  a  single, 
spatially  distributed,  strong  interfering  CW  signal. 
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The  research  results  summarized  in  section  V  above  have 
demonstrated  that  potentially  large  processing  gains  may  accrue 
to  a  narrow  band  HF  antenna  array  system  which  employs  adaptive 
beamforming  techniques  to  eliminate  a  single,  CW  interfering 
signal.  Adaptive  processing  methods  similar  to  those  used  in 
sonar  and  seismic  systems  have  direct  application  to  such  HF 
systems.  No  problems  unique  to  the  HF  environment  were  found 
to  influence  the  adaptive  processing.  The  performance  of  the 
HF  adaptive  beamforming  system  was  accurately  predictable 
using  theoretical  results  obtained  in  the  acoustical  signal 
processing  area. 

This  research  also  demonstrated  that  continuous  adaptive 
beamforming  is  desirable.  It  was  shown  that  "freezing"  the 
adaptive  process  for  as  little  as  60  seconds  can  cause  a 
degradation  of  8  to  10  db  in  output  signal- to-noise  ratio. 

Time  constant  investigations  indicated  that  as  few  as  680 
adaptions  are  required  to  achieve  a  3  db  improvement  in  output 
SNR  for  an  8  element  array  with  40  adaptive  weights.  The  al¬ 
gorithm  employed  to  obtain  these  results  requires  one  multiply 
and  add  for  each  weight  adapted.  Thus,  a  total  of  27,200  such 
operations  are  required  for  each  3  db  improvement.  Modern  low 
cost  digital  hardware  can  easily  perform  one  multiply  and  add 
in  two  usee.  The  time  constant  of  adaptation  would  then  become 
54.4  milliseconds,  a  rate  sufficiently  fast  to  track  ionospheric 
variations . 

In  addition,  this  research  program  also  suggested  and  in¬ 
vestigated  a  technique  for  implementing  stable  adaptive  digital 
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feedback  filters.  No  other  methods  are  presently  available  for 
feedback  filters.  Althouqh  only  preliminary  results  have  been 
obtained,  they  are  encouraging  and  may  lead  to  the  development 
of  an  adaptive  processina  system  for  use  in  arrays  with  a  large 
number  of  elements  and/or  simultaneous  beams. 


vi  r . 


sugcpstions  for  futurf  rusk arch 


The  research  results  and  conclusions  described  in  this  report 
have  demonstrated  that  adaptive  beamfomir.a  methods  have  great 
potential  in  HF  antenna  arrays.  Unfortunately,  these  results  are 
based  upon  a  data  sample  which  is  not  representative  of  the  system 
parameters  used  in  most  HP  applications.  To  this  end,  further 
research  should  be  conducted  to  examine  the  potential  of  adaptive 
methods  in  an  environment  which  bettor  approximates  that  of  actual 
systems.  The  specific  recommendations  are  as  follows: 

1.  To  record  and  process  data  observed  simultaneously 
at  array  elements  using  receiver  bandwidths  greater 
than  1  KHz. 

2.  To  investigate  the  performance  of  adaptive  beamforming 
methods  in  multiple,  hiqh  amplitude  interference 
environments  such  as  that  generated  by  teletype  trans¬ 
mitters  . 

3.  To  measure  the  capability  of  an  adaptive  processor  to 
track  time-varying  interference  --  i.e.  a  noise  source 
with  time-varying  frequency  spectrum  and/or  direction  of 
arri va 1 . 

4.  To  determine  the  effectiveness  of  adaptive  processing 
in  the  presence  of  siqnal-oonerated  clutter  and/or 
multipath.  Previous  work  has  been  based  on  the  assumption 
that  the  time  waveform  of  the  interference  is  uncorrelated 
with  that  of  the  desired  signal. 

r>  •  To  develop  adapt  ve  algorithms  better  suited  to  those 
problems  specific  to  HF  systems  --  c.u.  multiple  beam 
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processors,  ranqe/dopplcr  adaptation,  and  direction 
findimi  methods. 

6.  To  investigate  the  hardware  requirements  for  the  imple¬ 
mentation  of  real  time  adaptive  processors.  This  study 
should  include  an  investigation  of  the  relative  merits 
of  an  all  digital  processor  vs.  a  hybrid  analog/digital 
system. 

7.  To  further  explore  and  develop  procedures  for  utilizing 
feedback  in  adaptive  filters.  This  work  should  include 
array  studies  as  well  as  one-dimensional  filters. 
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APPENDIX 


THE  design  of  adaptive  tapped- delay  line  processors 

WHICH  EMPLOY  FEEDBACK 

Non-recursive  adaptive  tapped- delay- line  processors  of  the 
type  shown  in  Fig.  22  do  not  employ  feedback  signal  paths.  As 
a  result,  the  impulse  response  of  such  filters  is  necessarily 
Limited  to  the  total  length  of  the  delay  line.  In  order  to 
obtain  fine  frequency  resolution  with  this  processor,  it  may  be 
necessary  to  use  an  inordinate  number  of  taps  in  the  line. 

Since  the  convergence  rate  of  adaptive  processors  decrease  as 
the  number  of  taps  increases [22 ] ,  a  multiplicity  of  taps  may 
impose  prohibitively  long  convergence  times  on  the  adaptive 
processor;  particularly  in  the  case  of  array  processing  where 
N  lines  are  required  for  an  N  element  array. 

An  alternative  tapped-delay- line  processor,  which  does 
•  -rnploy  feedback,  is  shown  in  Fig.  45.  This  processor  is  commonly 
referred  to  as  a  recursive  digital  filter.  Because  of  the 
presence  of  feedback,  an  infinite  length  impulse  response  can 
be  obtained  using  a  finite  number  of  filter  coefficients. 
Alternatively,  a  recursive  filter  can  provide  an  equivalently 
fin-'  frequency  resolution  with  fewer  filter  coefficients  than 
for  a  non-recursive  design. 

One  problem  which  occurs  with  recursive  filter  design  is 
tha*-  of  obtaining  a  stable  impulse  response.  It  is  well  known 
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that  when  the  poles  of  the  feedback  polynomial  lie  outside  the 
unit  circle* ,  the  resulting  filter  is  unstable.  Since  the  feed¬ 
back  coefficients  determine  the  location  of  these  poles  in  a  non 
linear  manner,  it  is  not  an  easy  task  to  reflect  the  stability 
requirement  in  terms  of  the  coefficient  values.  For  this  reason 
simp]  *  adaptive  feedback  processors  have  not  yet  been  proposed. 

The  approach  taken  in  this  research  was  to  design  a 
minimum-phase  non-recursive  digital  filter  and  then  to  use  this 
filter  to  set  the  feedback  coefficients  in  a  recursive  design. 
Minimum  phase  filters  are  constrained  to  have  all  zeros  in  the 
o table  region  and  thus  the  resulting  recursive  design  necessari¬ 
ly  lias  all  poles  in  the  stable  region  and  is  therefore  stable. 

S i mo ] o  adaptive  procedures  are  available  for  the  design  of 
i  l.umum  phase  filters  and  have  been  widely  developed  in  the 
field  of  inverse  or  deconvolution  filtering. 

A.  DECONVOLUTION  AND  INVERSE  FILTERING 

A  filter  {g^}  is  defined  as  the  inverse  of  a  second  filter 
{hk  1  if  it  "undoes"  the  effect  of  {hk } ,  that  is,  for  any  input 
sequence  {W^ } 

[{W,.}*{l>k)l*{9k)  -  {Wk}  . 

"  f  *■'  is  often  referred  to  as  a  deconvoiution  filter 

★ 

The  -’.-transform  used  herein  is  defined  as 

X(?)  =  T  x  zK 
k  K 

using  a  positive  power  of  z  convention.  The  region  of  causal 
stability  is  z  >  1. 
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tor  (hk }  •  Takiny  z- transforms  of  both  sides  with 


this  means 


G(z)  =  z  {yk }  and  H(z)  =  z  {h.  } 


H(z)G(z)  --  1. 


G(z)  =  l./H(z) 


Such  a  require!  *nt  appears  simple,  but  has  several 
inherent  problems.  Most  serious  is  the  problem  of  stability. 

U  thG  fiitur  transform  H(z)  were  to  have  a  zero  inside  the 
unit  circle,  that  is,  if  it  were  non-minimum  phase,  the  required 
GU)  would  h  ave  a  pole  inside  the  unit  circle.  The  resulting 
{V  would  be  causal  but  unstable.  Secondly,  a  finite  duration 
impulse  response  fh^  would  invariably  require  an  infinite 
Jurat  ion  Sk )  as  its  inverse.  This  indicates  that  {g  }  could 
no  alwajo  be  realized  by  a  tapped  delay  line  or  non-recursive 
filter.  For  ease  of  implementation  such  a  non-recursive  or  all 
zero,  realization  is  desirable.  Thus,  subjected  to  such  con¬ 
straints,  there  may  exi^t  no  ideal  inverse  filter  {g  } 

To  avoid  such  problems  it  becomes  necessary  to  loosen  the 
definition  given  by  Eq.  (Al) .  Instead,  to  be  an  inverse,  it 
is  only  necessary  that  'gk}  do  the  best  possible  job  of 
inverting  hk }  .  The  performance  measure  to  optimize  is  the 
mean-square  error;  the  solution  that  minimizes  the  mean-square 
•  r  i  >1  is  the  Wiener  fi  lter.  In  its  non-recursive  form  it  is 
u  'fined  by  this  familiar  system  of  normal  equations 
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G  (N  )  p 

9.  '  RXX  -xd 


(A2) 


where 

g(N)t=  tg0gig2- - '  the  N  filter  coef£ioients  {gk>' 

P1'  .  [E(dvxk  .)]  0  <  j  <  N-l,  cross  correlation 

xd  K  K  J  t 

of  the  filter  input  )  and  the  desired  output  {dR)  ,  and 

Rxx  =  tE(xk-l  *k-jn  "i1, 

the  matrix  of  input  < xR >  autocorrelation.  For  the  special 
case  of  inverting  a  filter  <hR)  ,  where  {dR)  =  (WR)  , 

Pxd  *  lh0  0  0  •••  01 


and 


R. 


XX 


=  [E(hk_.  0  <  i»  j  i  N_1 


D.  TOEPLITZ  MATRICES 

When  dealing  with  stationary  random  processes,  the 
coefficient  matrix  of  the  system  of  normal  equations  Rxx  has 
a  pleasing  property:  it  is  Toeplitz,  so  that  all  elements  along 
any  diagonal  are  equal,  and  it  is  symmetric.  It  is  therefore 
defined  by  the  elements  of  its  first  row, 

rxx<l>  •••  txx<N‘1)1 


XX 
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in  1947  outlined  a 


where  r  (j)  ~  E(x,  x,  .).  Levinson^6^ 
xx  k  k-j 

simple  algorithm  for  solving  a  system  of  linear  equations  having 

a  Toeplitz  coefficient  matrix.  The  algorithm  requires  only  the 

first  row  of  coefficients,  thus  reducing  the  required  storage 
2 

space  from  N  to  N ,  over  conventional  Gaussian  matrix 
solutions  techniques.  Likewise,  the  number  of  operations  re¬ 
duces  by  a  factor  of  -oughly  N. 

An  additional  advantage  is  that  Levinson's  algorithm  uses 
a  recursive  approach?  it  solves  a  system  of  N  equations  in 
N  unknowns  using  the  solution  to  a  subsystem  of  N-l  equations 
in  N-l  unknowns.  By  starting  with  the  trivial  case  N-l 

rxx(0)90  =  p0 

the  algorithm  can  recurse  through  all  values  of  N.  Note  that 
the  mean-square  error  performance  can  be  calculated  at  each  step 
of  the  recursion.  When  the  error  is  small  enough,  the  recursion 
process  will  stop. 

C.  STABILITY  OF  INVERSE  FILTERS 

It  is  interesting  to  note  that  there  are  no  restrictions 
on  the  original  filter  when  using  Eq.  (A2)  to  find  the  inverse 
filter.  Recall  that  the  first  definition  given  for  inverse 

G(z)  =  l./H(z) 

required  ^hut  li(z)  =  zfh^}  have  no  zeros  inside  the  unit  circle. 


90 


r  1 

Ilore,  however,  since  {gk  }  is  derived  non-recursively  and  has 
only  N  coefficients,  its  time  response  must  be  stable  for  any 
(h^ }  •  It  might  be  asked,  what  good  reason,  if  this  stability 
is  artificial,  due  onxy  to  the  finite  duration  of  the  impulse 
response;  as  N  grows  indefinitely  will  instability  result? 

It  can  be  shown  that  {gk  }  will  indeed  be  stable  if  unrestricted 
in  length,  irrespective  of  {hk }  .  Furthermore,  as  N  grows 
(gk)  becomes  a  better  inverse  of  {hk }  in  the  mean-square 
error  sense,  and  if  ll(z)  is  minimum  phase, 

lim  Gv  (z)  =  1/H (z)  , 

N>n, 


which  is  the  intuitive  definition  set  forth  originally. 
If  H ( z )  is  minimum  phase  with  a  simple  zero 

H(z)  =  1  -  z/a  | a  j  >1 

This  result  indicates  that 


or  for  finite  N 


lim  G(N)  (z) 

N  -*oo 


r  (N )  ,  , 
G  (  Z , 


1_  _ 1 

hTz)  1-z/a 


1 

l-z/a 


Thus  the  N-l  zeros  of  the  non-recursive  inverse  (9^}  are 

distributed  as  to  approximate  a  pole  in  the  mean-square  error 

[27] 

sense.  Read  and  Burrus  1  showed  that  a  simple  pole  can  be 
approximated  by  N  zeros  distributed  uniformly  around  a  circle 


of  radius  a  with  the  zero  at  z  =  a  removed.  This  approxi¬ 
mation  was  derived  simply  from  the  truncated  impulse  response 
of  a  filter  having  a  simple  pole, 


1 

1-z/a 


z  2 

<!>  + 


+ 


N-l 


N 


Such  an  approximation,  however,  will  not  have  minimum  mean- 
square  error.  The  filter  G^(z)  derived  non-recursively  as 
above  has  its  zeros  slightly  perturbed  from  the  circle  z  =  a. 
For  example,  consider  the  case  of  a  pole  at  z  *  2  and  N  -  6 
shown  in  Table  Al. 

TABLE  Al.  INVERSE  APPROXIMATIONS 


TRUNCATION 

INVERSE 

OF 

(1  - 

( z/2 ) ) 

IMPULSE 

ZERO 

IMPULSE 

ZERO 

RESPONSE 

LOCATIONS 

RESI  7NSE 

LOCATIONS 

1.00000 

-2.00 

.99982 

-2.1719 

. 50000 

2 . 00e^ 12  0 

.49954 

2 . 1444e^121, 7 

.25000 

2.00e“^120 

.  24904 

2.1444e"j121,7 

.  12500 

2 .  00e3  60 

.12305 

2.0666e^62*3 

.03125 

2 . 00e~^  60 

.05860 

2.0666e"j62,3 

.02344 

Since 

(N ) 

G  ' (z)  approximates  1/H(z) 

in 

the 

mean-square 

error  sense 

,  it  follows  that 

1/G(N) (z) 

must 

similarly  approxi- 

mate  H(z).  This  leads  perhaps  to  the  most  important  property 


92 


■ 


f  If  I 

of  the  inverse  filter  proved  by  Me«jol*  .  For  .1  polynomial 

c(t,|U)  .  z;,k) 

the  ix>tj  must  always  lie  outside  the  unit  circle  in  the  complex 
z  plane,  if  (g^ )  is  the  minimum  mean-square  error  inverse  of 
any  filter  (h^  )  .  This  means  that  the  poles  of  a  filter  defined 
as 

H(N)(x)  -  l./G(fi)(z) 

must  be  outside  of  the  unit  circle.  With  the  poles  outside,  the 
new  filter  II  (z)  must  have  a  stab  It*  causal  time  response. 

D.  INVERSE  FILTERS  IN  RECURSIVE  FILTER  DESIGN 

The  non-r^cursivc  filter  is  most  popular  for  data  processing 
applications.  The  greatest  advantage  is  its  ease  of  design; 
while  it  is  extremely  difficult  to  optimize  a  recursive  filter 
in  all  but  the  most  structured  problems,  the  non-recursive  filter 
is  uniquely  defined  by  a  set  of  linear  equations.  Also,  it  ha3 
guaranteed  stability.  Since  the  non-rocursive  filter  by  defini¬ 
tion  has  no  feedback,  its  transfer  function  can  have  no  poles, 
only  zeros  in  the  plane,  and  thus  must  be  stable.  Furthermore, 
by  allowing  the  filter’s  delay  line  to  lengthen ,  its  impulse 
response  will  approach  that  of  the  optimal  causal  filter,  and 
its  moan-square  error  will  monotonically  converge  to  the 
optimal  performance.  Thus,  the  optimal  causal  filter  which  is 
recursive  may  be  approximated  to  any  degree  of  accuracy  by  the 
optimal  non-recursive  filter,  given  enough  delay  elements. 
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A  sonous  disadvantage  of  the  non-recursi vo  filter  configu¬ 
ration  is  its  poor  bandpass  capabilities.  For  estimating  a 
narrow-band  signal  in  noise,  a  good  filter  must  be  selective 
with  a  sharply  defined  passband.  With  no  poles  in  the  non- 
recursive  filter's  transfer  function,  a  narrow  passband  can  only 
be  approximated  by  placement  of  a  finite  number  of  zeros.  A 
recursive  filter  with  a  single  pole  pair,  on  the  other  hand, 

Will  in  general  have  a  better  bandpass  characteristic  tt  ei  a 
low  order  non-recursi ve  filter.  A  desirable  filter  would 
possess  the  selectivity  of  a  recursive  filter,  while  retaining 
the  design  advantages  of  the  non-recursive  filter. 

The  transfer  function  of  a  recursive  filter  takes  the 


form 


H  ( z) 


A  (2) 

qTzT 


where  A (z)  and  3(z)  are  polynomials  in  positive  powers  of 
Such  a  filter  will  be  unstable  if  any  of  its  poles,  that  is, 
zeros  of  B ( z) ,  fall  inside  the  unit  circle  in  the  z  plane. 

However,  if  both  numerator  and  denominator  were  inverted  non- 

(N)  (N) 

recursively  with  C  (z)  the  inverse  of  B(z),  and  D  (z) 


the  inverse  of  A(z),  the  ratio 

r,(N) 


II 


(z)  = 


C(N>(z> 


D(N>(z) 


z . 


would  necessarily  be  stable,  since  D(z)  has  its  zeros  outside 
the  unit  circle  according  to  Szego's  theorem.  Therefore,  to 
design  a  recursive  filter  by  this  technique,  it  is  only 
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necessary  to  find  the  expression 


,  .  _  A(z) 

H(z)  -  but 


such  that  H ( z)  approximates  an  optimal  filter  Hn(z)  and  such 


that 


r.(N)  _  C(N)  (z) 

11  {t) - fNl 

(z) 


approximates  il(z)  sufficiently  well.  Note  that  if  H(z) 
were  both  stable  and  minimum  phase,  so  that 

lim  C(N)(z)  =  1/B  ( z )  and  lim  D(N)(z)  =  1/A(z) 


then  it  follows 


,  .•  ..  (N)  1/B  (  z )  _  A  (  Z )  _ 

^  "  lz)  I7aTTT  -  bTIT  -  H(z)  ’ 

i: .  DBS  IGN  APPROACHES  AND  EXAMPLES 

Throe  different  design  approaches  were  taken.  The  first 
two  used  the  fact  that  the  unconstrained  optimal  Wiener  filter 


takes  the  form 


4>  (z) 


H  .  (z)  =  ~ 
opt  <t> 


for  the  estimation  problem,  with  the  signal  power  spectrum 


4SS(Z)  =  J  «3S(k)  2 

k  =  -™ 
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and  the  input  power  spectrum 


*xx(z)  *  J  *xx<k)  z 

k=-  oo 


l^ach  infinite  summation  is  truncated  and  then  used  to  form  an 
approximation  to  the  optimal  filter 


M 


Hj(z)  = 


k=-M 


I  *ss(k)z 


2M 


I.*ss(k"M)2 


k=0 


M  . 

I  ^xx(k)z 

k=-M  xx 


2M 

l  *xJk-M)^ 

k»0  xx 


_  A(z) 
k  "  B<z>  ' 


Finally,  both  A(z)  and  B(z)  are  inverted  by  (z)  and 

(N) 

C  (z)  respectively  as  explained  earlier  to  form  the  stable 
filter 


h^n)  (Z)  =  )  (z) 

D(N)(z) 


(A3) 


Till s  is  called  the  double-sided  autocorrelation  inversion  filter. 
The  second  approach  was  similar,  except  that  the  truncation  took 
only  positive  powers  of  z,  giving 


HII (z) 


*  k 
I  ♦  s(k)zK 

k=0  SS 

M  v 

1  >vX(k)Z 

k  =  0  xx 


A  ( z) 
B(z) 


(A4) 


This 


led  to  the  single-sided  autocorrelation  inversion  filter. 
The  third  approach  is  considerably  different.  Let 
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*ss(z)  '  Ma(*)lla(i) 
*xxU)  -  "x(2,llx<j) 


With 


Hs(z) 


Both 

Then 


filters  II 

s 

the  optimal 


(a)  and  1^(2)  are  minimum  phase 
causal  filter  takes  the  form 


and 


stable. 


Hoptc 


V*> 

H75T 


Wl'oru  thB  additional  polynomial  0(21 
2ur°  l0catl0ns  defined  by  the  causality 
da  approximation  to  this  expression  be 


accounts  for 
constraint. 


any  new 
Next  let 


H<2) 


H  (z) 

i-yir 


Wi  tli 

!,s(z) 


r(N)  . 

z)  inverting  »x(2)  and  D 
t*le  inversion  processor  is  then 


(N )  ,  . 

(z)  inverting 
given  by 


HIII (z) 


r(N) 

^TnT 


(z) 

(z) 


( A5) 


A1J  three  filters 
ir volving  a  narrowband 
overall  signal  to  noise 
he  seen  in  Figs.  46  and 


were  tested  using  a  specific  problem 
signal  omleddod  in  broadband  noise, 
ratio  was  -lldb.  The  power  spectra 
47.  In  addition  to  the  three  filter 


The 

can 
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0 

Fig. 


46. 


0,125  0.250  0.375 

Normalized  Frequency  (Hz) 

Power  Density  Spectrum  of  Desired 


0.500 

Signal 
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Power  Desnity  (db) 


,i"ai-"S  ',,-Sor,‘,-d  ‘,bov"-  . .  HI  tors  wo  ro  designed: 

tU"  U,,CO,>3trai"'  U  f'Hor.  Uv  causal  Wionur  filter,  and 

tho  non- recursive  Wiener  filter.  The  logarithmic  magnitude 

responses  of  each  filter  is  shown  in  Figs.  48-53.  Note  that 
each  of  the  five  recursive  filters  is  far  superior  in  its 
bandpass  selectivity  to  the  single  non-recursive  filter. 

The  im  an-square  error  evaluation  indicated  that  the 
-wly  derived  filters  were  inferior  with  the  error  defined  as 

ek  =  <SR  '  V2 

=  filter  output 

TO  evaluate  tho  performance  independent  of  time  delay  between 

actual  signal  value  and  the  estimated  value,  the  error  was 
redefined  as 


2  o 

e,  =  (s,  -  V  \ z 

k  k-D  k' 


Uiiu  B  was  allowed  to  verv  tHp  , 

y.  ine  universal  minimum  for  each 

filf^r  was  then  found  TaVnim  at 

Und*  Table  A2  summarizes  these  results. 


Table  A2.  MEAN-SQUARE  ERROR  CALCULATIONS 
TYPE  OF  FILTER  B  =  0  VARIABLE  D 

i)  WIENER  UNCONSTRAINED  .511  >5il 

‘  *  CAUSAL  .683  .611 

3)  NON-RECURSIVE  .  750  .734 

4)  DOUBLE-SIDED  AUTOCORRE-  .9R4  coc 

NATION  Eq.  (A3 )  *585 

5)  SINGLE-SIDED  AUTOCORRE-  .812  cQ1 

CATION  Eq.  (A4 )  *591 

b)  AVERSION  PROCESSION  92  1 

hq.  (A5 )  *632 


OUTPUT 

SNR 

-1.803  DB 
-4.384  DB 


0.614  DB 


too 


Amplitude  Response  (db) 


Fig.  48.  Amplitude  response  for  optimal  non-causal 
filter . 
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Amplitude  Response  (db) 


Normalised  Frequency  (db) 


Fig.  49.  Amplitude  response  for  optimal  causal 
filter. 


Amplitude  Response  (db) 


Amplitude  response  (db) 


1 


Fig.  51.  Amplitude  response  for  suboptimal  processor 
Hj  (Z)  given  in  Eq.  (A3)  . 
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Amplitude  Response  (db) 


For  narrowband  signals  in  broadband  noise,  the  recursive 
filters  derived  using  the  properties  of  inverse  filtering  may 
perform  quite  well,  if  dilay  between  the  signal  estimate  and 
the  actual  signal  value  is  not  important.  Indeed,  in  many 
applications  even  delays  of  several  seconds  may  not  be  signifi¬ 
cant.  The  technique  used  in  the  derivation  of  these  filters  is 
a  new  way  of  imposing  stability  on  a  recursive  filter  where 
before  instability  was  a  problem.  It  could  probably  have 
greatest  impact  in  an  extension  of  the  work  set  forth  here, 
the  problem  of  adaptive  recursive  filtering. 

F.  DISCUSSION  AND  CONCLUSIONS 

Encouraging  results  have  been  obtained  indicating  that 
minimum  phase  filters  can  be  used  successfully  in  the  design 
of  recursive  Wiener  filters.  Three  specific  design  procedures 
were  suggested  and  investigated.  Each  design  is  amenable  to 
the  implementation  of  adaptive  methods.  Further  research  into 
this  problem  in  continuing  to  determine  the  applicability  to 
array  processing. 
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